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Abstract 


Using an approach free from momentum extrapolation, we calculate 
the nucleon magnetic moment and the fraction of the nucleon spin carried 
by the quark angular momentum in the quenched lattice QCD approxi¬ 
mation. Quarks with three values of lattice masses, 210, 124 and 80 MeV, 
are formulated on the lattice using the standard Wilson approach. At ev¬ 
ery mass, 100 gluon configurations on 16 3 x 32 lattice with fd = 6.0 are 
used for statistical averaging. The results are compared with the previ¬ 
ous calculations with momentum extrapolation. The contribution of the 
disconnected diagrams is studied at the largest quark mass using noise 
theory technique. 


*This work is based on V.G.’s research for a dissertation to be submitted to the Graduate 
School, University of Maryland, by V.G. in partial fulfillment of the requirement for the PhD 
degree in physics. 
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I. INTRODUCTION 


The magnetic moment is one of the fundamental properties of the nucleon. Exper¬ 
imentally this quantity has been measured to a very high precision [11]. In the nucleon 
model-building, the magnetic moment is usually the first to check against the experi¬ 
mental data. To understand how well a particular version of lattice QCD can simulate 
the internal structure of the nucleon, the magnetic moment is an observable that one 
must study after the nucleon mass. 

Theoretical computations of the nucleon magnetic moment using lattice QCD tech¬ 
niques have been performed by a number of groups on lattices of different physical 
sizes and lattice spacings with different nucleon sources (interpolating fields) |2[- 0. 
Most calculations done so far have used the quenched approximation. And almost all 
calculations begin with the lattice calculations of Sach’s magnetic form-factor Gm(Q 2 ), 
where Q 2 = — + g 2 , at discrete values of the momentum transfer q. The magnetic 

moment is just Gm{Q 2 ) at Q 2 = 0 which can not be obtained more directly in this 
approach because the relevant off-forward matrix element is linearly proportional to 
q. The lattice magnetic momentum is obtained usually by extrapolating these discrete 
values at finite Q 2 to 0. 

However, there are potential problems in the Q 2 extrapolation, mainly due to the 
finite volume effects. On a finite lattice, the momentum transfer q = (q \, g 2 , gQ is 
quantized: g* = 2n-n/Ni in unit of inverse lattice spacing, where iV* is the total number 
of sites in momentum direction i and n = 0,1,..., iV* is an integer. The smallest possible 
non-zero Q 2 is equal to — q% + (27r/Max(IVj)) 2 , where g 0 = \J M 2 + (27r/Max(A^j)) 2 — 
M, and M is the nucleon mass. For a lattice with a spatial dimension 16 3 , Q 2 m in = 


2 M(yj (vr/8) 2 + M 2 —M). For (3 = 6, this corresponds to ~ 0.5 GeV 2 . At this scale, the 
nucleon form factors are very sensitive functions of Q 2 . Indeed Ge(Q 2 ) increases from 
Q 2 = 0.5 GeV 2 to Q 2 = 0 by almost a factor of 3. Therefore, an extrapolation across 
this range of Q 2 must be strongly model dependent if without any prior knowledge of 
the Q 2 dependence of the nucleon form factors [4|.pI1 |TTf| . Fortunately, the experimental 
data shows that the Q 2 -dependence of the nucleon form factors can be fitted to a 
dipole form. Moreover, in a limited Q 2 range, the dipole form is not too different from 
a monopole. Therefore, we find in the literature both forms are commonly used to 
fit the lattice data |^.|T^]TI]]. An example of the extrapolation using different fitting 


formula can be found in Ref. [pj]. 

In this paper, we report a more direct calculation of the nucleon magnetic moment 
using the elementary definition in electromagnetism [Q. Although this approach was 
first mentioned in Ref. [Q, to our knowledge, however, no systematic study along this 
line has ever been reported in the literature. As is quite obvious, the great advantage 
of this method is that it is free from the extrapolation of the finite momentum transfer 
Q 2 . We will comment on a different finite volume effect later in this paper. 

Using a similar approach, we can calculate the fraction of the nucleon spin carried 
by the quark angular momentum. In recent years, the study of the spin structure of 
the nucleon has stimulated much interest in both the experimental and theoretical nu¬ 
clear and particle community [|T3|,|n|]. It was found in Ref. | [Hj| | that the quark angular 
momentum in the nucleon can be obtained through deeply-virtual Compton scattering 
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and other hard exclusive processes. While this observation has stimulated much per¬ 
turbative QCD and phenomenological study [16], the quark angular momentum can 


also be calculated in lattice QCD. A first attempt in which the momentum extrapola¬ 
tion was used in the magnetic moment calculation was reported in Ref. [[L7j]. As shown 
in Ref. 


15 


it can be calculated using the direct theoretical definition as well. In this 
paper, we also report such a calculation. 

Our paper is organized as follows. In Sec. II, we give a brief theoretical discussion 
of the nucleon magnetic moment and quark angular momentum in the continuum. We 
state our conventions and describe the implementation of the continuum formulas in 
the lattice simulations. In Sec. Ill, we outline the setup of our lattice calculations and 
present our main results. We will compare them with previous calculations. We also 
consider the contributions of the disconnected diagrams to both the magnetic moment 
and quark total angular momentum. In the final section, we present the conclusion 
and discuss possible ways to improve upon the present calculation. 


II. THEORETICAL DETAILS 

In this section, we present a number of theoretical formulas which are useful in our 
numerical calculations. In the process, we will make clear our notations and conven¬ 
tions. 


A. Magnetic Moment and Quark Angular Momentum 


According to the standard electromagnetic theory the magnetic moment operator 
of a system can be defined as 0: 


h = 2 / ^ x .7em]d 3 r , 


( 1 ) 


where j em is the electromagnetic current density. Its matrix element in a quantum state 
\jrn z = j) (( jm\jm) = 1) quantized in the ^-direction defines the magnetic moment, 


h = 0 = j\frz\jm z = j) . 


( 2 ) 


We will use this definition to calculate the magnetic moment of the nucleon in its rest 
frame. 

The magnetic moment can also be obtained from the form factors of the electro¬ 


magnetic current. For the nucleon, we have |jj 

(P'S'lj^mPS) = u(P'S') \yFi(Q 2 ) + i 


2m n 


f 2 (Q 2 


u(PS) , 


( 3 ) 


where a tw = £[ 7 ^, 7 "], | PS) represents the ground state of the nucleon with four- 
momentum P and polarization S, q = P' — P is the momentum transfer, Q 2 = —q^ + cf, 
and u(PS) is an on-shell Dirac spinor. Fi(Q 2 ) and ^(Q 2 ) are the well-known Dirac 
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and Pauli form factors. The nucleon magnetic moment can be obtained from Sach’s 
magnetic form factor 


G m (Q 2 ) = F^Q 2 ) + F 2 (Q 2 


(4) 


at Q 2 = 0: fijy = Gm( 0). This approach has been used in most of the lattice QCD 
calculations. 

To understand the spin structure of the nucleon, one needs the QCD angular mo¬ 
mentum operator in a gauge-invariant form [ |T5| ,p~4[ 

<4qcd = Jg + Jq , (5) 


with 


J.= 


d 3 r[f x T] — / d 3 r 


^—'i/j + ^ [f x (— iD)]%b 


■J„ — I d 3 r 


f x [E x B\ 


( 6 ) 


where S is the Dirac spin matrix and D = V — is the covariant derivative, 
T = (T j0 ) is a component of the quark energy-momentum tensor to be defined later. 
In a nucleon state with hclicity 1/2 and momentum along the ^-direction, the total 
hclicity can be calculated as the expectation value of J z , 


- - -S + Lq + Jg - J q + Jg 


(7) 


where the various contributions to the spin of the nucleon are defined as the expectation 
values of the corresponding operators. For instance, the matrix element of J q defines 
the quark total angular momentum contribution to the nucleon spin. 

As in the case of the nucleon magnetic moment, the quark contribution to the nu¬ 
cleon spin J q can be obtained from the form factors of the quark energy-momentum 
<—> 

tensor T^ u = (the parentheses around the two indices means symmetriza- 

tion and subtraction of the trace and D=D — D) 


(P’S'\T^(0)\PS) = 
u(P'S‘ 


p'V'ritQ 2 ) + 1 t 2 (q 2 ) + At-U(q 2 ) 


2 M 


M 


u(PS) , 


( 8 ) 


where P M = (P M + P lf1 )/ 2. The quark contribution to the spin of the nucleon is then 
simple J q = |[Ti(0) + T 2 (0)]. In Ref. |l^ , this was used to calculate J q . 


B. Lattice Formalism 


Light (up and down) quarks are put on the lattice using the standard Wilson 
formulation. In the discrete Euclidean spacetime, the fermionic part of the QCD action 
is given by Hl9f : 
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FIG. 1. Two-point and three-point diagrams. 


S w . (9) 

where i,j is the sum over the full set of spacetime, Dirac, color indices with i = {x, a, a} 
and j = {y, /3, b}. M %3 is defined as 

M u = K X) {(! _ 7 n)ocyU^{x) ab 5 X:y _j X + (1 + r y tl ) a pUj l (x - y) ab 8 x ,y +fl } , (10) 


where k = 1/(2 (mo + 4)) with mo the bare quark mass, y is a unit vector along the 
H = 1,2, 3,4 direction, and the lattice spacing a is set to 1. U^x) = e lA ^ represents 
a link connecting the nearest lattice sites going from x to x + fi. Our notation for 
7 -matrices follows Refs. 


2C.21 


7f = 7° = 





( 11 ) 


where E means Euclidean space and will be omitted henceforth. The standard contin¬ 
uum Dirac field ^ is related to the above by an extra factor of yJm,Q + 4. 

We are interested in calculating two- and three-point correlation functions involving 
the nucleon interpolating fields on the lattice. The graphical representation for them 
is shown in Fig.|T|, where to is the “time” of the nucleon creation (source position), t x 
is the “time” of the nucleon annihilation (sink position) and t y is the “time” of the 
operator insertion (we sum only over the spatial coordinates x and y). 

The nucleon two-point correlation function is defined as 

G 2 aa'(t x ,P') = E e ^(°IXa(^(0)|0) , (12) 


where x = {t xi x) are discretized spacetime coordinates, a, a' are Dirac indices, and 
|0) represents QCD vacuum. Operator y is an interpolating held with the nucleon 
quantum number. We choose j|77| pi| 

Xa = e» k (u‘C^d‘)u k a . (13) 


Indices {i, j, k} = 1, 2, 3 represent colors, C = 7472 is a charge-conjugation matrix. The 
time dependence of the two-point function can be obtained by inserting a complete set 
of intermediate states | n) with the nucleon quantum number: 
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( 14 ) 


G 2aa '(t x ,P’) = E(°IXa(0)|n)(n|x^(0)|0)^e- E " ' 

n 

where we have used the Euclidean translation xiPPx) — e^ Htx ~ lx ^x{G)eG Htx+lx ^ and 
V 3 is the three volume. In the limit t x ^>> 1, only the lightest state (namely, the nucleon) 
with momentum P' contributes. Therefore the nucleon mass can be extracted from 
the (large) time dependence of the two-point function, 


M = hi (15) 

G 2aa (t x + 1, P' — 0) 

M is in the unit of the inverse lattice spacing. I 11 the next section, we will extract the 
nucleon mass this way as it is needed for the magnetic moment calculation. 

The matrix element of an observable O in the nucleon state can be extracted from 
the following three-point function: 

G M (t x ,t y ,P,P') =Y.^ iP ' 3 e >S (0\x«(x)dWx a ,{0m ■ (16) 

S,y 

At large time separation, to <S< t y <<C t x , it can be written as 
Gsaa' (hr, t y , P, P') = 

E 3 2 J2(°\Xa(0)\P'S')(P'S , \d(0)\PS)(PS\x a ,(0)\0) ■ e~ Ep,tx • e -(^e-®e')h ; (17) 

S,S' 


where we have kept only the intermediate nucleon state. The labels S and S' denote 
the nucleon polarization. The nucleon state here is assumed to be normalized to 1. 

I 11 this paper, we are interested in the electromagnetic current operator j^iy) 
and quark energy-momentum tensor operator T^ u (y). The continuum expression 
j e ™(y) = ei'0i7/W ; * (where ip is the quark held of flavor i and e* is the electric 

charge) is a conserved current . On a lattice, the above expression can be imple¬ 

mented straightforwardly, but it is not conserved because of the finite lattice spacing. 
Instead, is multiplicatively renormalized if the power-suppressed contributions are 
neglected. The conserved electromagnetic current in the Wilson fermion formalism has 
the following more complicated form |27] : 


= \ [V’O/Xtm- l )Uv.{y)P(y + A) + + m)(ta* + l ) u l(y)P(y)\ • ( 18 ) 

Differences between local and conserved currents on a lattice are discussed in Refs. 

The quark energy-momentum tensor T /U ,(y) is not conserved. Therefore, there is 
no preferred way to put the operator on a lattice. We choose the local form, T^(y) = 
|'0(y){7 M [D — D] v }ip(y) as in the continuum with the derivative D ^ defined on a 
lattice as: 


= 2 Pp{x)il>(x + p)~ Ul(x - p)ip(x - p ,)] 


(19) 
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The above definition is related to the continuum one, say, the MS scheme, through 
a finite renormalization. Ignoring the mixing contribution from the gluon energy- 
momentum tensor, the tadpole improved renormalization constant for operator T^ u 
has been calculated perturbatively and is near unity, Z = 1.045 [Bj. A nonperturbative 
renormalization technique may be needed to find the renormalization factor reliably. 

Finally, we consider the transformation of operators from the Minkowski to Eu¬ 
clidean space. When the vector current is defined with the Euclidean 7 matrices, 

the corresponding Lorentz vectors in Eq. (3), for example, must be defined with 
V 4 = V° and Vi = —iVfa. Alternatively, one can define the current as then the 

Euclidean four-vectors are related to the Minkowski ones by V 4 = iV° and Vi = V^. 


C. Formulas for Physical Observables on Lattice 


As we have discussed earlier, the nucleon magnetic moment and quark total angular 
momentum can be calculated by an extrapolation of the relevant form factors to Q 2 = 0 
0- H. On a lattice, these form factors can be extracted from the ratio of three- to two- 
point correlation functions. For the magnetic form factor Gm(Q 2 ) = Fi(Q 2 ) + F 2 (Q 2 ) 
and the total angular momentum T(Q 2 ) = \\Ti(Q 2 ) + T 2 (Q 2 )\, one can obtain at the 
large time separation, 


E €*»-*'*{ 0 \Xa(x)Ji(y)xAom 

x,y _ 

Ee-^(0|x 7 (x)x y (0)|0) 

X 


f 


(P r + M) 7r 


-e ijk a j 

aa 


,q k G M {Q 2 


( 20 ) 


Ee 

x,y 


iqy-iP'x^Q | (a .) T 4i (y)xA o) 10 > 

£e-^0|x 7 (Ew(0)|0> 


P° 


/ 


i (P' + M) 7r 


-e ijk a j 

aa 


q k TlQ 2 


( 21 ) 


where a is the Pauli matrix, and the multiplier / is equal to 


Epi T M 

t y 

tx 

2Ep 

l _y_ 

tx 

’Ee-^(0| Xa (x)x Q ,(0)|0)' 

X 

2Ep> 


_Ep T M_ 


E e _ '^(0 x«(s)x a /( 0 ) 0 ) 


X 


f = 1 when q = 0. These expressions are used in this study to calculate form factors 
to check against the existing results. 

The nucleon magnetic moment and quark total angular momentum can be calcu¬ 
lated by taking the derivative over the momentum transfer in the limit q —» 0 : 


T,(Q\xa(x)yjJi{y)x a ’(fy\Q) , Ee 

X,y _ _ ° X,y 

E(0|x 7 (z)x 7 '(0)|0) d(iqj) 


d f 

d(iqj ) (P f + M) 77 , 


iki k 
6 (7 

^ ^ aa‘ 


,q J G m (Q 2 ) 


2 \ Q 


iqy—iP'x 


(0|x«(a;)Ji(s/)x a '(0)|0) 


E e~ iP ' s (0 1 y 7 (x)x , ( 0 ) 10 ) 


*0 iV jk a k c 

2 M 


N , 


( 22 ) 
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(23) 


E^IxaMsjiMfo^tojio) Ef* ip '*(o|xt.Wr4i(s)x c ,'(o)|o) 

S,y _ _ O g,y _ 

£(0|x 7 (z)%y(0)|0) 9 (^i) Ee-^(0|x 7 (x)x 7 /(0)|0) 

X X 


d ifP° 
d(i qj ) (A + M) lY 


e ijk a 


k 

aa' 


T(Q 2 ) 


?-► o 


e ijk a 


k 

aa 
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The summations over x and y ensure the nucleon has vanishing three-momentum 
(in the rest frame) and the forward matrix elements are calculated. Therefore the Dirac 
index 7 ' must be the same as 7 . The nucleon can be polarized in the different spatial 
direction by selecting different a and a'. The new results obtained from the above 
formulas will be presented in the next section. 


III. NUMERICAL RESULTS 

I 11 this section, we present our numerical calculations on the 16 3 x 32 lattice at 
(3 = 6.0 with the Wilson formulation of fermions. The results include the nucleon 
mass, form factors, magnetic moment, and the quark total angular momentum. 

The coupling (/?) corresponds to lattice spacing a ~ 0.1 fm or a' 1 ~ (1.7—1.9) GeV. 
The antiperiodic boundary condition for the Dirac operator is used in the time di¬ 
rection. We use 100 quenched configurations at every quark mass to attain reason¬ 
able statistics. However, as the quark masses approach the chiral limit (k reaches 
K cr ~ 0.1568(1) [[H]]), the statistics require many more lattice configurations, and hence 
the error bars increase markedly for a fixed number of configurations. 

We use three different values of the mass parameter k = 0.152, 0.154, 0.155, cor¬ 
responding to quark masses 210, 124, 80 MeV, respectively. Point sources are used 
to create the nucleon: interpolating field \ a i s placed on a well-defined initial lattice 
point, called 0, and twelve different sources with different color and Dirac indices are 
used to start the conjugation gradient. 

Errors are determined by the standard jackknife procedure [^8|- [pi]]. We average the 
three-point and two-point functions separately using this method for 100 configurations 
and put already averaged values in the final formulas for the magnetic moment and 
quark total angular momentum. 


A. Nucleon Mass 


As is clear from the formulas in the last section, we need the lattice nucleon mass 
to calculate the nucleon electromagnetic form factors and magnetic moment. We use 
Eq[L5| to extract the nucleon mass on the lattice. The result is shown in Fig^j as a 
function of the time-slice for k = 0.152. A stable plateau is seen to have reached 
between the time slices 10 and 17. 

Table | shows the comparison between the nucleon masses obtained in our calcu¬ 
lations and those in the literature (Ref. PUTTH). The results, shown in the units of the 
inverse lattice spacing a -1 , are in good agreement. Note that the lattice nucleon masses 
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FIG. 2. Nucleon mass as a function of time-slice, calculated on 16 3 x 32 lattice with 100 
quenched gauge configurations at f3 = 6.0 and k = 0.152. 


are consistently bigger than the physical mass of ~ 0.534 in the lattice unit because of 
the large lattice quark masses. 


B. Nucleon Magnetic Moment and Quark Total Angular Momentum 

We first present the results for the nucleon form-factors at several momentum trans¬ 
fers Q 2 and compare them with the previous calculations []5J. The expected consistency 
is a useful check of the codes for three-point functions. Results are compared in Tables 
H] (for q 2 = 0.154) and flj (for q 2 = 0.308) and are seen in agreement within error bars. 
The agreement improves for larger ft. One noticeable trend is that when the quark 
mass becomes smaller, the magnetic form factor decreases. This could be an indication 
that the finite lattice size effect is important. We will comment on this further below. 

Results for the nucleon magnetic moment are summarized in Table m One obvious 
conclusion from the table is that the lattice magnetic moments are about 30% smaller 
than the experimental data (/i p = 2.79 and = —1.91). There are a number of factors 
which can account for this discrepancy. One is the use of the quenched approximation. 
We will report a full dynamical calculation in a future publication. Another is the 
lattice quark masses. The smallest quark mass in our calculations is 80 MeV for 
ft = 0.155, which corresponds to the pion mass ~ 0.5 GeV. In a phenomenological 
study 0,^], Leinweber et. al. demonstrated that as the pion mass approaches zero, 
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TABLE I. Nucleon mass in lattice units ( ma ), calculated on 16 3 x 32 lattice with 100 


quenched gauge configurations at [3 = 6.0 


K 

nucleon mass 

nucleon mass from Ref. [4 

nucleon mass from Ref. [ 171 

0.152 

0.86(2) 

0.87(3) 

0.882(12) 

0.154 

0.7(3) 

0.73(5) 

0.738(16) 

0.155 

0.63(5) 


0.67(15) 


TABLE II. Magnetic form-factors Gm(q 2 ) f° r (fa 2 = 0.154, calculated on 16 3 x 32 lattice 
with 100 quenched gauge configurations at j3 = 6.0 


K 

proton 

neutron 

proton from Ref. ||] 

neutron from Ref. [£i|] 

0.152 

1.36(9) 

-0.88(7) 

1.22(7) 

-0.781(59) 

0.154 

1.13(12) 

-0.738(75) 

1.17(11) 

-0.748(69) 

0.155 

1.02(15) 

-0.661(101) 




TABLE III. Same as Table [0], but for fa 2 = 0.308, calculated on 16 3 x 32 lattice with 
100 quenched gauge configurations at f3 = 6.0 


K 

proton 

neutron 

proton from Ref. ||] 

neutron from Ref. [[|] 

0.152 

1.09(8) 

-0.697(65) 

0.906(59) 

-0.586(47) 

0.154 

0.877(97) 

-0.567(93) 

0.895(90) 

-0.578(89) 


TABLE IV. Nucleon magnetic moments, calculated on 16 3 x 32 lattice with 100 quenched 


gauge configurations at (3 = 6.0 


K 

proton 

neutron 

0.152 

2.24(19) 

-1.46(13) 

0.154 

2.09(30) 

-1.38(21) 

0.155 

1.88(37) 

-1.26(28) 
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TABLE V. Quark total angular momentum, calculated on 16 3 x 32 lattice with 100 


quenched gauge configurations at (5 = 6.0 


K 

connected J u +d 

disconnected Ju+d 

0.152 

0.47(7) 

-0.12(6) 

0.154 

0.45(9) 


0.155 

0.44(10) 



the magnetic moment shows a steep rise as a result of chiral singular contributions. 
In our calculations, we observe an opposite trend: the magnetic moment decreases 
along with the lattice quark masses. This is seen in the calculation of the form factors 
at finite momentum transfers and a number of other observables |35|]. We suspect 


strongly that this is a finite volume effect which will go away when we increase the 
physical dimension of the lattice by 50% or so. 

The results for the quark total angular momentum contribution to the nucleon spin 
are shown in Table [V]. In this case, our results are quite consistent with those in Ref. 
17| , although we have avoided the use of the extrapolation in Q 2 . For instance, our 


result J u+ d = 0.44(10) at k = 0.155 from the connected diagram can be compared with 
J u+ d = 0.44(7) at k ci = 0.1568 in [[L7]]. The leading chiral contribution to this quantity 
was studied in Ref. [^3| and is not strong. As the quark mass approaches the chiral 
limit, the lattice result shows a slight decrease, although the large error bars prevent a 
clear-cut comparison. 

Although the quark angular momentum from the connected diagram accounts for 
nearly all of the spin of the nucleon, there are other contributions which cannot be 
neglected. For instance, the sea contribution through the disconnected diagrams is 
expected to quench the connected result, resulting a total quark angular momentum 
accounting for about half of the nucleon spin ||34| . 


C. Disconnected Diagrams 


The contribution of the disconnected diagrams represents one of the possible correc¬ 
tions to the nucleon magnetic moment and quark total angular momentum presented 
in the previous subsection. The disconnected insertion is shown in Fig^j. Disconnected 
diagrams are difficult to calculate in general because they involve the closed-loop con¬ 
tributions resulting from the self contractions of the quark fields in an observable. An 
efficient method to obtain the matrix elements is the noise theory discussed in Refs. 
36|- 


39fl , explained as follows. 


Consider a set of random sources rj l (x) satisfying 
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FIG. 3. Disconnected diagram. 


1 ^total 

—- r] ] \x)r]\y) = S xy . (24) 

1 ’total i = i 

The total number of sources IVtotai is a key parameter controlling the desired statistics. 
In this work, the random sources rf(x) are created by putting 1, — 1, % or —i randomly 
on every lattice point for every Dirac and color index. Three-point functions for the 
disconnected diagrams can now be expressed as: 

1 -^total 

GfTAt„t y ,P') = G M (t x ,P') ■ — Y. Yri V (y)[yxd(y)]M-\y,z)r,\z) , (25) 

-('total i =l 


where O = j cons or T j0 , G 2 aa '(t x , P') is Green’s two-point function, and M^ x (y,x) is a 
Wilson-fcrmion propagator. 

Unfortunately, this method requires a large N totSL \ to obtain reasonable statistics, 
which is very expensive in practical simulations. It was discussed in Refs. p8| , |39| that 
the so-called expansion method can expedite the calculation significantly. The idea is 
to approximate the full fermion propagator, 


M (x, y) = 


$xy KjP xy 


(26) 


where P xy = Y^A+-y y W y A)5x,y-y + (l--y y )UUx-fj,)S xx+fI } by one that includes only 
the first few terms in the small n expansion, 


^exp(u y) — $xy + K-Pxy + K? P xy + ••• • (27) 

One then writes M~ l (x,y) = M e “p((c, y) + {M^(x, y) — M^(x, y)}. The first term, 
M^(x, y), can be calculated directly without conjugation gradient, and the remainder, 
M l {x, y) —M~ x p(x, y), is computed by the noise method discussed above. Because the 
remaining contribution is small, the number of random sources needed for the required 
statistics is considerably reduced. 

In this study, we keep ten terms in M e “The remainder is calculated with 30 
random sources. The results are obtained only for the smallest hi = 0.152 because even 
in this case the disconnected contribution to the magnetic moment has almost 100% 
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error. For the quark total angular momentum, the result is more stable and the error 
is about 50%. Our final result for the disconnected diagram is J u+d = —0.12(6). This 
agrees with [[17]] J u+d = —0.094(24). The negative contribution cancels part of the 
result from the connected diagram. We have no simple explanation for the sign. 


IV. CONCLUSION 


In this paper, we report a lattice QCD calculation of the nucleon magnetic moment 
and the fraction of the nucleon spin carried in the quark total angular momentum. We 
use the standard Wilson formulation for fermions and 100 quenched configurations in 
the Monte Carlo evaluations of the Feynman path integrals. 

Our results for the magnetic moment cannot be compared directly with the ex¬ 
periment yet. We are uncertain, for example, about the systematic error caused by 
the quenched approximation, large lattice quark masses, and the finite volume effect. 
The smallest quark mass in our calculations is about 80 MeV (k = 0.155) which is 
about ten times larger than the physical quark mass. The pion mass dependence of 
the magnetic moment has been studied recently in Refs. |fTT].[32|| . In these papers, a 


strong dependence on the pion (quark) mass in the region m^ 2 < 0.5 GeV is suggested. 
A phenomenological chiral extrapolation of the previous lattice results to the physical 
pion mass gives results consistent with experimental values. 

One interesting finding of this study is that the magnetic form factors, magnetic 
moment, and the quark total angular momentum decrease with the quark masses. The 
most plausible explanation is that when the quark mass becomes lighter, the nucleon 
size grows, and the present lattice volume cannot fully accommodate a physical nucleon. 

We have also studied the contribution of the disconnected diagrams to various 
physical observables. The error bars grow very significant at small quark masses. The 
sea contribution to the quark total angular momentum is found to be negative. The 
fraction of the nucleon spin carried by the up and down quark angular momentum 
is J q = 0.35 A 10 at s = 0.152. The result for a fully dynamical simulation will be 
reported elsewhere. 
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